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Type I and Two-Gap Superconductivity in Neutron Star 
Magnetism 
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ABSTRACT 

Neutron-star inner cores with several charged baryonic components are likely to be 
analogues of the two-gap superconductor which is of current interest in condensed- 
matter physics. Consequently, type I superconductivity is less probable than type II 
but may nevertheless be present in some intervals of matter density. The intermediate- 
state structure formed at finite magnetic flux densities after the superconducting tran- 
sitions is subject to buoyancy, frictional and neutron-vortex interaction forces. These 
are estimated and it is shown that the most important frictional force is that produced 
by the stable stratification of neutron-star matter, the irreversible process being diffu- 
sion in the normal, finite magnetic- flux density, parts of the structure. The length-scale 
of the structure, in directions perpendicular to the local magnetic field is of crucial 
importance. For small scales, the flux comoves with the neutron vortices, as do the 
proton vortices of a type II superconductor. But for much larger length-scales, flux 
movement tends to that expected for normal charged Fermi systems. 
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1 INTRODUCTION 

The evolution of neutron star magnetic fields has been a 
topic of considerable interest since the early papers of Pacini 
(1967) and Gold (1968), the discovery of pulsars by Hewish 
et al (1968), and the paper of Baym, Pethick & Pines (1969) 
on proton type II superconductivity. There has been much 
uncertainty in the magnetic flux transport properties of 
those parts of the interior where the density does not ex- 
ceed the nuclear density po = 2.5 x 10 14 g cm -3 . For a 
typical neutron star of mass IAMq, the central density pre- 
dicted by many equations of state is p ~ 2 — 4po, indicating 
the presence of hyperons and of even greater uncertainties 
in magnetic flux transport properties in the inner core. It is 
evident that empirical deduction of these properties is not 
feasible and that a priori theoretical input is required. 

The general assumption has been that the protons form 
a So type II superconductor. As functions of matter density, 
calculated proton energy gaps are typically A p ~ 0.5 MeV at 
p ~ po but decrease rapidly at higher densities (see Pethick 
& Ravenhall 1995; Heiselberg & Hjorth- Jensen 2000). Thus 
it seems likely that the condition for type II superconductiv- 
ity, k — A/£ > 1/a/2, where A is the penetration depth and 
£ the coherence length, will be satisfied at p ~ po, that is, in 
regions of the outer core where matter is composed of pro- 
tons and x So or 3 P2 superfluid neutrons, charge-neutralized 
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and in weak-interaction equilibrium with electrons and neg- 
ative p-mesons. Calculations of A p are even less reliable at 
higher densities, p ~ 2po, so that the possibility of type I su- 
perconductivity in the inner core cannot be excluded. There 
has been some work on E~-hyperon pairing (see Vidana & 
Tolos 2004) which indicates a large So gap and the possi- 
bility that the inner-core superconductor has two or more 
components. But much uncertainty remains concerning all 
baryonic gap calculations, and some authors have questioned 
the existence of type II superconductivity at any core den- 
sity (Link 2003; Buckley, Metlitski & Zhitnitsky 2004; but 
see also Alford, Good & Reddy 2005). 

A further possibility is that unconfined quarks are 
present and that the long-term stability of the core mag- 
netic field is then determined by the properties of the Meiss- 
ner effect for a colour superconductor (Alford, Berges & 
Rajagopal 2000). If this very interesting question is to be 
studied by observations on neutron star magnetic fields it is 
essential that there are analyses of magnetic flux transport 
in more prosaic systems, such as baryonic type I and two- 
gap superconductors. The present paper is addressed to this 
problem. 

Magnetic flux transport for type II superconductivity, 
assuming the matter composition anticipated at p ~ po, 
has been analysed previously (see Jones 2006; also refer- 
ences cited therein). The problem is well-defined owing to 
the existence of the mixed state in which, on microscopic 
scales, magnetic flux is quantized in units of <j>o = hc/2e = 
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2.07 x 1(T 7 G cm and confined to the cores of proton vor- 
tices. Type I superconductors behave differently at the tran- 
sition from the normal state in the presence of a magnetic 
field owing to their positive surface energy, equivalent to the 
condition k < l/y/2. The intermediate state formed mini- 
mizes surface area as much as external constraints allow. 
It consists of flux-free superconductor in equilibrium with a 
filamentary structure of normal protons and magnetic flux 
density B w H c , where H c is the thermodynamic critical 
field. Averaged over volumes large compared with the scale 
of the structure, the magnetic flux density (B) is that of the 
initial normal system. Owing to the external constraints, the 
minimization condition does not, in itself, lead to structures 
that are well-defined or of universal form. Being of little 
practical utility, type I superconductors have been less well 
studied than type II. We refer to Tinkham (1996) for further 
details. 

As in the type II case, we assume that an approxima- 
tion to static hydrodynamic equilibrium exists in the neu- 
tron star interior at times before the superconducting phase 
transitions. This equilibrium changes because the spatially- 
averaged components of the type I superconductor stress 
tensor are larger than those of the normal-state Maxwell 
tensor by a factor of H C /{B), which is of the same order as 
for type II superconductors (Jones 1975, Easson & Pethick 
1977). The change in the stress tensor occurs in a time 
short compared with any possible flux expulsion time be- 
cause cooling of the neutron star interior is rapid at the 
superconductor critical temperature. The problem is to de- 
termine how the magnetic flux filaments move under any 
buoyancy force which may then appear. 

The presence of S~ hyperons at p ~ 2po would make 
the neutron star interior an example of the two-gap su- 
perconductor which is of current interest in studies of liq- 
uid metallic hydrogen at high pressures (Babaev, Sudb0 
& Ashcroft 2004). Given that several baryonic components 
may be present in the inner core, it might be thought that 
its superconducting properties are complex. But it appears 
that the complete system can be a type I superconductor 
only if all individual Fermi liquids are either normal or type 
I. 

Both the two-gap case and the factors that determine 
the form and scale of the filamentary type I structure are ad- 
dressed in Section 2. The various forces acting on a moving 
type I filament are found in Section 3. That arising from the 
stable stratification of neutron star matter (see Goldreich 
& Reisenegger 1992, who introduced the concept as a factor 
limiting non-solenoidal ambipolar diffusion of magnetic flux) 
is by far the most important in all parts of the core where 
the negative /i-meson threshold is exceeded. In these regions, 
filament motion is limited by the rate of lepton diffusion be- 
tween the filament and the flux-free superconductor. The 
order of magnitude of the force between a neutron vortex 
and a type I filament is evaluated in Section 4. It is shown 
that neutron vortices are effectively trapped by the struc- 
ture, as is the case in type II superconductivity. Section 5 
gives a summary of the circumstances under which the inner 
core could retain magnetic flux for times long compared with 
the ages of most pulsars or the spin-down times of the neu- 
tron stars in binary systems. If these can be excluded, any 
evidence of long-term flux freezing would indicate the pres- 
ence of more exotic core structures. In common with the 



previous paper on type II superconductivity (Jones 2006) 
the scope of this paper is intended to be purely technical 
and it does not consider observational evidence relevant to 
magnetic flux evolution. 



2 TYPE I AND TWO-GAP 
SUPERCONDUCTORS 

The intermediate state of a type I superconductor has no 
universal form (see Tinkham 1996) and, in a neutron star 
core, must be largely defined by the flux distribution in ex- 
ternal regions where it can be regarded as frozen for the 
short time-scales relevant to the superconducting transi- 
tion. Given that constraint, we anticipate that the spatially- 
averaged magnetic flux density (B) will initially be identical 
with the normal-state flux distribution existing before the 
transition. We assume that, in the initial stages of the tran- 
sition, small flux-free regions form with a growth rate gov- 
erned by ohmic diffusion in the surrounding normal phase. A 
cylindrical shape would minimize surface area and least per- 
turb the surrounding magnetic flux distribution. Unless (B) 
is quite close to H c , growth of these flux- free regions must 
lead to their coalescence and the formation of a filamentary 
structure of magnetic flux and normal protons occupying a 
fraction (B) / H c of spatial volume. 

The temporal order of the various superconductor tran- 
sitions which occur is uncertain. An isothermal neutron star 
core cannot be assumed at the early times in which these 
transitions occur because its composition may permit direct- 
Urea processes in some regions but not others. However, it 
is possible to infer that the filamentary structure is unlikely 
to be a two-dimensional lattice. That it is probably inter- 
connected in a complex way can be seen from the following 
considerations. Suppose that, at any instant, type I super- 
conductivity is confined between spherical surfaces of radii 
n,2 (t) with T2 > r\. As the neutron star cools, these sur- 
faces move because they are dependent on the evolving tem- 
perature distribution T(r,t) and on the critical temperature 
T c (r, t). Filaments with a distribution of cross-sectional area 
and therefore varying total flux form between these surfaces. 
The case that ri < must produce complexity because it 
involves the merging of two independently formed filament 
distributions subject to flux conservation. (In the simple case 
of a pure dipole field, the distributions merging would be 
those formed on the two sides of the magnetic equatorial 
plane.) There would be complete complexity for n(oo) = 
and partial complexity for n(oo) > 0. This time-dependence 
of n would be consistent with the baryon-number depen- 
dence of calculated proton energy gaps at inner-core den- 
sities (see Pethick & Ravenhall 1995; Heiselberg & Hjorth- 
Jensen 2000). Complex interconnections are also formed as 
a result of the fusion of adjacent filaments which move, un- 
der the buoyancy force, with velocities that depend on their 
radii (see Section 3). Our assumption about filament orien- 
tation is that the spatially-averaged magnetic flux density 
(B) will initially conform with the pre-transition flux dis- 
tribution. Any large-scale rearrangement at the time of the 
transition would be energetically not allowed. 

The filaments are subject to a post-transition buoyancy 
force, fs per unit volume. Its origin is more simple than in 
the type II case owing to the macroscopic cross-sectional 
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area of the filaments. The pre-transition Maxwell tensor 
expressed in terms of the spatially-averaged fields contains 
products (Bi)(Bj). These are replaced by {BiBj} in which 
B w H c . During motion toward a new equilibrium, the cross- 
sectional area of a filament may change owing to flux conser- 
vation and p-dependence of A p and there may be inward or 
outward lepton or proton diffusion, constrained by electrical 
neutrality. 

If superfluid £~ hyperons are present (the threshold is 
at p ~ 2po in many equations of state), there is negligi- 
ble Josephson coupling with the protons and the system is 
analogous with the two-gap superconductor of current in- 
terest, the particular example being that of liquid metal- 
lic hydrogen at high pressure (see Babaev et al 2004). We 
can assume that the internal temperature satisfies T <C T c 
throughout most of the time interval in which the flux dis- 
tribution evolves. Thus the complicated phase diagram of 
such systems does not concern us and the behaviour is sim- 
ple. The reasons for this will be given here in outline only 
because they follow very closely the work of Babaev (2002) 
on the structure of composite vortices, to which we refer for 
further details. For a system with several charged Fermi liq- 
uids i = 1,2..., the Ginzburg-Landau free-energy functional 
is written down in terms of the condensate amplitudes 
bare masses mi and charges e, of the components, and varia- 
tion with respect to the vector potential A gives the current 
density J. In the approximation that the |$»| are position- 
independent, the condition J = at radii s » A from a 
vortex gives a simple expression for A in terms of the gra- 
dients of the condensate phases on circular paths about its 
axis. (The total changes of phase x» on these paths are the 
phase windings.) The line integral of A gives the enclosed 
flux which, for two components with e\ — —e^ and opposite 
phase windings xi = — X2 of magnitude |xi,a| = 27r, is ex- 
actly 0o . This result, a single flux quantum, holds for any 
number of components provided e 2 and ei\% are both con- 
stants independent of i. Substitution of A into the Ginzburg- 
Landau functional shows that, for these phase windings, the 
long-range kinetic terms vanish identically and that single 
flux-quantum composite vortices are the lowest-lying states 
of the superconductor that can support magnetic flux. (The 
case of a zero phase winding for one of two components gives 
the interesting class of vortices with a fractional quantum 
of flux, but substitution of A into the free-energy functional 
shows that the long-range kinetic terms are then finite and 
that the vortex self-energy is much greater than for the sin- 
gle flux-quantum case. These vortices are not of interest at 
the temperatures T <C T c considered here.) It follows that 
a system of several components, at least one of which satis- 
fies the surface energy condition k > l/\/2, will behave as 
a type II superconductor even at magnetic field strengths 
that exceed the thermodynamic critical fields of all but that 
one component (see Babaev 2002). 

The conclusion is that if several charged baryonic com- 
ponents are present, as may be the case in the inner core, 
the system is a type II superconductor unless no component 
satisfies the k > 1/V% condition. This simple result indi- 
cates that type I superconductivity may be limited to no 
more than a small fraction of the core volume. 



3 THE FORCES ON A MOVING TYPE I 
FILAMENT 

3.1 Leptonic frictional forces 

Under a buoyancy force fs, the filamentary structure of a 
type I superconductor moves to some extent as a single en- 
tity because the nature of the stress tensor inhibits differen- 
tial motion and the resulting filament curvature. To estimate 
the forces that balance fs, we assume that the filaments 
are locally cylindrical, so minimizing surface area at con- 
stant flux, and are described by cylindrical polar coordinates 
s = (s, 9) moving with the filament. Except where otherwise 
stated, the radius s = a is assumed to be macroscopic, large 
compared with lepton gyro-radii inside filaments or with 
scattering mean free paths. For brevity, we consider bary- 
onic matter composed of protons and superfluid neutrons 
only. Then the relevant mean free path is that for electron 
scattering by muons. We consider, initially, frictional forces 
that arise from interaction between filament and leptons. 

Easson & Pethick (1979) have given the transport relax- 
ation time for electrons charge-neutralized by normal pro- 
tons, number densities iV e = N„. It is 
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In this expression, a is the fine structure constant; ftFe and 
€f p are the electron Fermi wave number and proton Fermi 
energy. The screening wave number is determined by the 
proton Fermi wave number k Fp - For nonrelativistic protons, 
it is, 



k F r 



( Aak Fp m p c\ 



(2) 



The appropriate adaptation of equation (I) for scattering 
by nonrelativistic muons, in the presence of protons, with 
independent number densities N e / iV M =fc N p is, 
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This is valid for Np such that k F ^ k F T, and the opti- 
mum definition of kFT is that given by the most massive 
particles in the system, the superconducting protons. This 
relaxation time is long (of the order of 10~ 12 s at T — 10 s 
K), so that the mean free path is always many orders of 
magnitude greater than the electron gyro-radius which is 
r B = € Fe /eB = 3.3 x lO^B^ 1 cm, for e Fe = 100 MeV, 
where B12 is the magnetic flux density in units of 10 12 
G. The appropriate transport relaxation time t p for elec- 
tron scattering by protons in the presence of nonrelativistic 
muons is given by an expression identical with equation (3) 
except that proton kinematic variables replace those for the 
muon. 

The assumption made here is that the boundary condi- 
tion satisfied by the lepton fluid velocity V; is (v;)x =0 at 
s = a on the filament surface. This differs from the condi- 
tion v;(o) = for conventional viscous flow which assumes, 
effectively, that a moving particle completely transfers its 
parallel momentum component to a surface on collision. 
This is the basis, for example, of the treatment of flow along 
a pipe under the Knudsen condition of kinetic theory (see 
Kennard 1938). However, there is negligible transfer of the 
parallel component in the case of lepton interaction with 
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a filament, effectively a cylinder of magnetic flux density 
B w H c . The small transfer which does occur is a conse- 
quence of scattering by leptons or protons. But the lepton 
mean free path is many orders of magnitude larger than the 
gyro-radius and becomes infinite in the T — > limit. On 
the basis of these considerations, we shall ignore the moving 
surface as a source of vorticity, which it would be in conven- 
tional viscous flow, and assume that the lepton flux Nivi 
is an irrotational and solenoidal vector. Its potential sat- 
isfies Laplace's equation and is determined by a Neumann 
boundary condition (see Batchelor 1967) on a static surface 
in coordinates fixed in the rotating star but instantaneously 
coincident with those of the moving filament. The velocity 
Vj is easily obtained for the case of constant number density 
Ni. In this frame, it has the familiar dipole form, with com- 
ponents (vi) s — —Ua 2 cos8/s 2 and (vi)e = —Ua 2 sm9/s 2 , 
for filament velocity U perpendicular to its axis. By inte- 
gration of the expression for the dissipation rate per unit 
volume given in terms of the stress tensor (see Landau & 
Lifshitz 1959), we find that the electronic viscous force per 
unit length acting on a filament is, 



f e = — 87r?7 e U, 

in which the shear viscosity is, 



(4) 



(5) 



(Easson & Pethick 1979) and is of the order of 10 20 g cm -1 
s _1 at 10 8 K. It is worth noting that the same calculation, 
made for a sphere of radius a, gives a force —12nri e a\J which 
is twice the Stokes force, an example of the minimum dis- 
sipation theorem (see Batchelor 1967, p. 227). Equation (4) 
is independent of a provided a ^> cr^ and so can be sig- 
nificant for thin filaments. For smaller values of a, but such 
that a^$> rs, the appropriate expression for the force can be 
found by an elementary classical kinetic theory calculation. 
The filament is regarded simply as a cylinder of magnetic 
flux density B w H c . The electron momentum and flux vec- 
tors are transformed into the filament rest frame and the 
momentum transfer obtained directly by integration over 
its surface. Under this condition, the force is given by, 
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(6) 



and is temperature-independent. Forces generated by inter- 
action with muons are given by expressions of the same form 
as equations (4)-(6) subject to replacement of electron num- 
ber density with that for the muons, kp e by Af^, and £p e 
by the muon kinetic energy £f m - For nonrelativistic muons, 
the relaxation time, 
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(7) 



replaces that in equation (5). 

Equations (4) and (6) have been obtained under the as- 
sumption of constant Ni but it is worth mentioning that a 
further frictional force arises if the scale length C for varia- 
tion of a leptonic number density, probably iV M , with depth 
is small compared with the neutron star radius. The rea- 
son is that the distributions of the fluid velocities v e and 
produced by filament movement are then not exactly iden- 
tical. Relative motion of the two fluids gives a force per unit 



length of filament of the order of, 
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It is negligible at small a compared with f e , and at large a is 
less significant than the force arising from the stable strat- 
ification of neutron star matter whose order of magnitude 
will be calculated here in much more detail. 



3.2 Stratification and stability 

It is assumed here that, immediately following the supercon- 
ducting transition and before filament movement, the con- 
ditions for charge neutrality and weak-interaction equilib- 
rium in terms of number densities and chemical potentials, 
N p — N e + iV M and fj, n — Hp = He = Hn + m M , are satisfied 
in the filaments and in the flux-free superconductor. With 
this definition of muon chemical potential, which excludes 
rest mass, the pressure can be expressed as a sum of its 
components, 

2 



P = Ph + Pe + P„ = Ph + ~N ef I e 



(9) 



where P h is the baryonic component and P e ,n are the lep- 
tonic components. We consider the effect of a small but finite 
outward displacement r — > r+8r in the position of a filament 
on the equilibrium of the matter in its interior. The initial 
problem is to determine the effect of the movement on the 
leptons and superconductor outside the filament. As noted 
in Section 3.1, the resulting particle fluxes Ar p>e>M v Pje , M , de- 
fined by a Neumann boundary condition at the surface of 
the moving filament, are irrotational and solenoidal. Thus 
the number densities N p ^ e ^ remain unchanged as functions 
of position in a frame of reference fixed to the rotating star 
even though weak-interaction transition rates are negligi- 
bly small. A displacement over the same distance 8r within 
the proton superconductor outside the filament is associated 
with number density changes N n — > N n (l — ("„) etc in each 
particle type and with a pressure change, 



8P = 



-PC - -PC 



(10) 



in which the coefficients are defined as Q n , P = 
Nn,p(dP/dN n ,p). Our assumption is that the equilibrium 
conditions are always satisfied everywhere in the supercon- 
ductor outside the filaments. Thus we can eliminate Cp.e.n 
and express SP as a linear function of £„, 



8P = 



Cn ^Qn ~\~ ^QpCpe + — Pe ~\~ —P[iC fie^J Cen^j , (H) 



in which the parameters are defined as follows; 

Me 

2^m' 



C, 



Cpe — 



N e + iV M C M e 



Cen — 



where, 



3 Me ApCp e 



d(nn - Hp) 



(12) 



(13) 



(14) 



(15) 



Neutron Star Magnetism 5 



and 



<9(/in - Hp) 

' dN v ' 



(16) 



But inside the displaced filament, the constraint of electrical 
neutrality, in the absence of diffusion and of significant weak- 
interaction transition rates, requires £ p = C, e = C,^ = £. (A 
tilde denotes quantities whose values inside a filament may 
differ from external values. The precise relation existing be- 
tween external and internal values at r would be replicated, 
under weak-interaction equilibrium, at r + Sr and so does 
not enter the problem. The effect of variation in H c over the 
displacement length Sr is too small to be significant.) The 
further constraints of pressure and neutron equilibrium are 
expressed as 8P = 8P and 5(i n = Sft n , respectively. Because 
fin is a function of both N n and N p in general, it is necessary 
to define a distinct internal value 7^ to represent the 
change in internal neutron number density. 

The pressure equilibrium constraint is then, 

8P = SP = -C„Q„ - C (Qp + \Pe + , (17) 

and the chemical potential constraint is, 

(^)^(C--C.) + (^)^ (Cp -c>o. (is) 

Equations (11), (17) and (18) can be solved for £ = G c Cn an d 
= GnCn- The differences between the chemical potentials 
inside the filament and those in the flux-free superconductor 
at the same radius r + 5r from the centre of the star are then, 



8/J-e — —/J-e (Gc — Cen) Cn, 



<W = jMm (G c 



Sup = (^) iV »( G »- 1 )C» 

dfip 



(19) 



(20) 



(21) 



The pressure increment <5P given by equation (11) can also, 
of course, be expressed simply in terms of the local matter 
density and gravitational acceleration g, 



8P — g (N n + N p ) m n 8r, 



(22) 



so that, from equations (11) and (19)-(22), the energy excess 
per unit volume of matter inside a filament can be expressed 
in terms of a force constant K, 



(Sr) 2 



(23) 



This equation defines K in terms of the composition- 
dependent parameters contained in (11) and (19)-(22), but 
the expression so obtained is not compact. Its evaluation is 
clearly very dependent on the choice of equation of state. 
The order of magnitude, for non-interacting Fermi gases 
(chemical potentials fi n — 128.2, fip = 8.2, fie = 120, pi M = 
14.3 MeV) gives K = 6 x 10 20 g 2 4 erg cm -5 . This is an inter- 
esting result. It shows that filament movement under a buoy- 
ancy force /s would be stopped in a distance Sr « fs /K ~ 1 
cm. 



This metastable state formed inside the filament can, in 
principle, convert to the equilibrium state by the irreversible 
loss of neutrino and thermal energy. But for superfluid and 
superconducting systems, baryonic semi-leptonic transition 
rates are negligibly small at T ~ 10 s K, unlike the nor- 
mal Fermi systems in the paper of Goldreich & Reisenegger, 
and they will not be considered further here. Direct leptonic 
transitions fj, ^ e with neutrino emission are excluded by 
the large difference between electron and muon Fermi mo- 
menta. The most significant process affecting the interior of 
a filament appears to be the diffusion of leptons and protons 
to or from the proton superconductor. These rates determine 
the filament velocity. 



3.3 Diffusion and filament velocity 

This diffusion, which was neglected in Section 3.2, means 
that the chemical potential differences given by equations 
(19)-(21) are actually functions of position, 5(j,i{s) for < 
s < a. The boundary condition Sfii(a) — is a consequence, 
for protons, of the strong interaction and, for leptons, of the 
fact that mean free paths in the flux-free superconductor at 
s > a are always many orders of magnitude larger than their 
gyro-radii at s < a. We are unaware of any full calculations 
of this diffusion process valid for B w H c and so have been 
obliged to rely on an order of magnitude estimate. We be- 
gin by noting that the field is marginally non-quantizing, as 
defined by Potekhin (1999), at T = 10 8 K and B = 10 14 G. 
A lepton typically completes many classical revolutions be- 
tween scatters and the gyro-radius is many orders of magni- 
tude greater than the Fermi wavelength. It then follows that 
the total electron and muon collision rates per unit volume 
are independent of B and are given by equation (3) and 
by an equation of the same form in which the muon kine- 
matic variables are replaced by those for the protons. We 
shall consider the electron case and assume that a scatter 
causes a classical guiding centre displacement of the order of 
rskFr/kFe, where tb is the gyro-radius. The radial electron 
flux is of the order of, 



r B k F T dSfi e \ ( r B k F T \ 

kFeknT ds J \ kFe ) 



(24) 



-D 



ds 



and defines the diffusion coefficient D, 

2 



D = 



He 

3k B T 



Tp Tf, j \ k F e / 



(25) 



which is a linear function of T. Diffusion is, of course, con- 
strained by the need to maintain electrical neutrality inside 
the filament. To allow for this, equation (24) should be mod- 
ified to include the radial electric field necessary to give iden- 
tical fluxes for leptons and protons, but we have neglected 
this problem and use the unconstrained electron diffusion 
rate given by equation (24) for the order of magnitude ve- 
locity estimate made here. 

The equilibrium considered at the end of Section 3.2 
continues to exist in the presence of diffusion because the 
chemical potential differences 5fj,i are maintained by out- 
ward movement of the filament. For a cylinder of radius 
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s < a inside the filament, the continuity equation is, 
f ^2ns'ds' = 2nsD (^) + ns'S, 



(26) 



in which the source term S depends on filament velocity U, 

(27) 



SB 



The condition for a time-independent <5/i e is, 
9<5/u, e _ sS 



(28) 



9s 2D' 

which, with the boundary condition 5p e (a) = 0, defines a 
steady-state velocity, 

4Df B 



U = 



Ko? 



(29) 



Evaluation for the parameters assumed in Section 3.2 gives 
a diffusion coefficient D = 2.1 x 10 _1 T 8 (H cli y 2 cm 2 s~\ 
where Ts is the temperature in units of 10 8 K and H c i4 is the 
thermodynamic critical field in units of 10 14 G. The velocity 

is, 



U — —gi4T$fB2t) (H C 14) 2 , 



(30) 



as we have argued in Section 3.1, is almost exactly irrota- 
tional. That it is not exactly irrotational and our neglect 
of quasi-particle excitations within the filament both imply 
that this cancellation cannot be exact. But these problems 
are not of significance for the results obtained in this paper. 



4 INTERACTION BETWEEN NEUTRON 
VORTICES AND TYPE I FILAMENTS 

It is generally assumed that, in regions of type II proton 
superconductivity and neutron superfluidity, the movement 
of magnetic flux is constrained by the strong and electro- 
magnetic interaction between neutron and proton vortices 
(Sauls 1989). Flux movement is also highly constrained in 
type I regions. This can be seen by calculating the change 
in free energy per unit length arising from the movement 
of a neutron vortex from the flux-free region of a type I 
proton superconductor to a position co-axial with a normal 
proton filament of magnetic flux density B w H c . We shall 
assume that the filament radius a is at least of the order 
of the neutron intervortex spacing (10 -3 — 1CP 2 cm). The 
free energy density for a system of neutron and proton su- 
perfluids can be expressed in the Ginzburg-Landau form as 



in units of cm s _1 . It is of a significant magnitude for filar a functional of the condensate amplitudes *„, p . It includes 



mentary structures with size less than a ~ 10 cm. 



3.4 The Magnus force 

All the forces considered above are frictional in character, 
and it remains to note the presence of a Magnus force on the 
moving filament. The proton superfluid circulates around a 
stationary filament of radius a at a velocity, 



v p = A. 

m„c 



(31) 



a term representing the phenomenon of superfluid entrain- 
ment, originally introduced by Andreev & Bashkin (1976) 
in connection with solutions of He 3 in He 4 , and further ex- 
tended to the magnetohydrodynamics of neutron-star super- 
fluids by Vardanian & Sedrakian (1981). It is likely that the 
inner-core neutrons form a 3 P2 superfluid, but for the or- 
der of magnitude evaluation made here we shall neglect this 
complication and assume 1 5o structure. In the approxima- 
tion that the \^ n , P \ are position-independent, it reduces to 
the intuitive form, 



We consider the case s>A, where A is the proton pene- 
tration depth, for which the vector potential can be chosen 



/£P i -f n i PP 2 . pn 

GL = fc + fc + Tjft V P + ft V P 



as, 



Ag = — H c Xexp 



. 1 nn 2 
Vn + -p 3 V n 

+ i-(VxA) 2 



(35) 



(32) 



For a filament moving with velocity U, the circulating super- 
fluid velocity is changed by the addition of a small increment 
<5v p . Flow is irrotational and solenoidal and is subject to the 
boundary condition v p ± — at s — a, so that the increment 
is 5v p g = 2U sin 6 at s = a in the filament rest frame. Ne- 
glecting entrainment terms, which do not contribute here, 
the isotropic component of the superfluid stress tensor, 



(see, for example, Alpar, Langer & Sauls 1984), in which /J'" 
are the proton and neutron condensation energy densities 
The superfluid velocities are, 

H 



% „ e 
75 V Xp A, 

2m v m„c 



(36) 



PP 

Ps VpiVpj 



(33) 



for superfluid density p p a p , is changed by a small increment 
—p pp v p eSv p e. Integration of this over the filament surface 
gives a transverse force of magnitude, 



/m = p7u 



(2TvaeH c X 
m p c 



(34) 



per unit length of filament in which, as expected, the final 
bracketed term is the circulation. 

This force is not large compared with /s for macro- 
scopic values of a and in any case, is cancelled by the lepton 
flow which is subject to the same boundary condition and, 



where \™,p are the condensate phases and m„ iP the bare 
masses. We refer to Alpar et al for the definitions of the 
superfluid densities p™"' pp and of the entrainment term p v s n 
in terms of the bare and effective particle masses. Variation 
with respect to A gives the electromagnetic current density, 

J = — (P? p v p + P f"Vn) • (37) 
m p 

It will be convenient here to re-use the cylindrical polar co- 
ordinates (s, 6) to describe neutron vortex structure. The 
condition J = at s >• A fixes the proton superfluid veloc- 
ity v p and, for zero proton winding phase, the asymptotic 
vector potential, 



p pn he m p 1 
p pp 2e m n s 



= Aq\/s 



(38) 
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whose line integral gives the fractional quantum of magnetic 
flux 4> described by Alpar et al. 

Movement of the neutron vortex from the flux-free re- 
gion of superconductor transfers its fractional flux quantum 
to the filament. There is a filament volume change which 
depends on the relative orientation of the vortex flux and 
the magnetic flux B in the filament. This depends both on 
the spin direction of the star and on the sign of the entrain- 
ment density pf n , and so is unknown. To be specific, we shall 
assume here, initially, the parallel case. Because the conden- 
sation energy density and H c are related by — = Hc/8tt, 
the filament free energy is increased by 4>H c /4n per unit 
length. There is no proton supercurrent inside the filament. 
Thus the total increase in free energy per unit length is the 
sum of this term and the change in the volume integral of 

fGL, 



AE = 



<j>H c 
4tt 



J 27VS ds (±p p /vl + pTv p ■ v„ + ± (V x A) 2 ) (39) 

in which the upper limit of integration is defined by the 
spacing between neutron vortices. Its order of magnitude 
can be evaluated for a model vortex defined by, 



Aos/X, 
A \/s, 



< s < A, 
6 > s > A, 
< s < £„, 
b > s > £„, 



(40) 



where £ n is the neutron coherence length and v® — 
h/(2m n (; n ). In this model, the vortex has flux (f> — \y\<f>o, 
where y = pJ n //°S P - Entrainment reduces vortex self-energy 
independently of the sign of y. For typical values (H c = 10 15 
G, In b/X — 10, £ n <C A) and the proton number density as- 
sumed in Section 3.2, the energy difference per unit length 
of vortex is AE = 1.6 x 10 7 |y| +4.6 x 10 7 j/ 2 erg cm -1 . The 
first term would become negative if the vortex flux and B 
were antiparallel. The parameter y is probably smaller than 
unity but is not well known, particularly at the high baryon 
densities of the inner core, so that AE could, in principle, 
be of either sign. But its magnitude is large, of the order of 
10 7 erg cm -1 or 0.6 MeV fm _1 . At this number density, the 
penetration depth is A = 8 x 10~ 12 cm and the interaction 
force, AE/X ~ 10 18 dyne cm -1 , is some orders of magni- 
tude larger than the Magnus force ~ 10 12 (v„) dyne cm -1 
arising from any plausible value of the spatially-averaged 
neutron superfluid velocity (v n ) relative to pinned vortices. 
Except for the possibility of sliding motion, flux movement 
is constrained in type I regions, as in type II. 

The interaction energy estimated here leads to con- 
clusions differing from those of Sedrakian, Sedrakian & 
Zharkov (1997) who have previously considered the equi- 
librium state of neutron vortices in a proton type I super- 
conductor with finite magnetic flux, parallel with the angu- 
lar velocity vector of the star. These authors assumed that 
the vortices are surrounded by coaxial normal filaments of 
radius a and uniform magnetic flux density, with flux-free 
superconductor at radii a < s < b, where b is determined by 
the intervortex spacing. They obtained a as a function of b 
and H c . However, entrainment reduces vortex self-energy so 
that the case AE > is the more probable. The stable equi- 



librium is then a state in which the vortices are surrounded 
by flux-free (except for the vortex core) superconductor at 
< s < a and the normal protons and magnetic flux arc 
confined between coaxial cylindrical surfaces, a < s < b. 



5 CONCLUSIONS 

In Section 2 we have observed, following recent work on 
two-gap superconductors (Babaev 2002; Babaev, Sudb0 & 
Ashcroft 2004), that a type I distribution of magnetic flux 
in a system of several charged baryonic components is pos- 
sible only if no superconducting component satisfies the 
k > 1/V2 condition. Type I characteristics are perhaps less 
probable than the ordered vortex-lattice type II structure 
for this reason. Nevertheless, it is possible that type I re- 
gions are present and the purpose of this paper is to see 
how magnetic flux distributions evolve within them. 

Our assumptions about the nature of the superconduct- 
ing transition have been stated in Section 2. A filamentary 
structure of magnetic flux and normal baryons forms oc- 
cupying a fraction (B) / H c of spatial volume. For a small 
window of possible (B) values, this fraction so large that 
the structure is inverted, the filaments being flux-free. In 
this case, the magnetic flux evolution would be that of a 
normal Fermi system and we refer to Goldreich & Reiseneg- 
ger (1992) for discussions of ambipolar diffusion and of the 
Hall effect. In particular, the force derived in Section 3.2 
is present above the muon threshold and removes the pos- 
sibility of ambipolar diffusion because the lepton diffusion 
considered in Section 3.3 is very slow over distances of the 
order of the neutron star radius. 

The movement of magnetic flux in a type I region is 
not, of course, independent of its behaviour elsewhere in the 
star (and vice- versa). This inconvenient mutual dependence 
is simply the result of flux conservation and of the large in- 
crease in energy that would be the consequence of indepen- 
dent movement in different regions. The problem is there- 
fore extremely untidy: estimates of the filament velocity U 
under the buoyancy force fs are made under the unknown 
constraints imposed by external regions. 

We consider the case in which (B) / H c is perhaps an 
order of magnitude smaller than unity. If the filament size 
distribution has very small values of a, it might be thought 
that the system should resemble an irregular type II vortex 
lattice with varying (large) integral numbers of flux quanta. 
We refer to Jones (2006) and references therein for further 
details of magnetic flux movement in a type II supercon- 
ductor. But there are some essential differences. It is not a 
lattice owing to the complex interconnections described in 
Section 2, and filament fusion is favoured because it reduces 
surface energy. The significant forces acting at small a are 
those given by equations (4) and (6). Here, rapid diffusion 
of leptons and protons reduces chemical potential differences 
to negligible values so that stable stratification has no effect. 
The buoyancy force per unit length of filament is 7ra 2 fs and 
it follows that the values of U are a-dependent. These dif- 
ferential values of U lead to contact, fusion of filaments to 
larger a, and increased values of U. But eventually, larger 
a and longer diffusion times cause stable stratification to 
become important. At an internal temperature T w 10 9 K 
appropriate for a young neutron star, the velocity U is deter- 
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mined by equation (30) for filament sizes larger than a ~ 1 
cm. Therafter, the progress of fusion continues throughout 
the type I region until it is stopped either by the external 
factors we have discussed above or by filament velocities be- 
coming too small. 

External factors are the more important constraint lim- 
iting the size of a. For macroscopic values, the decrease in 
surface energy through fusion is quite negligible compared 
with the increase in energy arising from any large-scale cur- 
vature or distortion of an individual filament. For example, 
under the buoyancy force fs alone, the equilibrium central 
displacement of a filament of length L, fixed at its ends by 
external constraint, would be limited and of the order of 
8nL fs / H c . This means that if neutron stars have type I 
and II regions in concentric radial shells, flux movements 
within them are not independent. 

Type I filamentary structures interact strongly with 
neutron vortices, most of the interaction energy AE aris- 
ing from superfluid entrainment. For small a, the structures 
move easily under the buoyancy force fs or a force fv derived 
from interaction with neutron vortices, as do proton vortices 
in the type II case. Reference to equation (30) shows that 
comovement with vortices at radial velocities of 10~ 7 — 10~ 6 
cm s _1 occurs easily. Velocities of these orders of magnitude 
are possible during the spin-down of young isolated pulsars, 
such as the Crab, or in the propeller- phase spin-down of 
neutron stars in binary systems. It is, however, unfortunate 
that our conclusion depends crucially on the value of a and 
on the diffusion coefficient D, of which equation (25) is no 
more than an order of magnitude estimate. For values larger 
than a ~ 10 3 - 10 4 cm, the frictional force produced by dif- 
fusion becomes large enough to limit these velocities. As a 
approaches its limiting value (~ 5 x 10 5 cm) the speed of 
flux movement tends to that expected in a normal proton 
system as described by Goldreich & Reisenegger (1992) and 
referred to earlier in this Section. 
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